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Abstract 

In this paper, we present a new hierarchy which includes the intermediate long wave (ILW) 
equation at the lowest order. This hierarchy is thought to be a novel reduction of the 1st 
modified KP type hierarchy. The framework of our investigation is Sato theory. 

1 Introduction 

It is well known that Sato theory was established by M.Sato around 1980 to give a unified 
viewpoint for integrable soliton equations [1]. A lot of important studies have been based on this 
magnificent theory since then. In this paper, we focus on one of the basic ideas of the theory, 
summarized as follows [2]. 

"Start from an ordinary differential equation and suppose that the solutions satisfy certain 
dispersion relations. Then, as conditions the coefficients must satisfy, we obtain a set of nonlinear 
partial differential equations. If we assume a particular set of linear dispersion relations, we obtain 
the KP hierarchy as the corresponding PDEs." 

It is known that most nonlinear partial differential equations which have ./V-soliton solutions 
correspond to certain equations of KP hierarchy. But, in the case of the intermediate long wave 
(ILW) equation, although its ./V-soliton solutions and an inverse scattering transform were well- 
known [3,4,5], the correspondence to the KP hierarchy has remained unclear. 

The ILW equation was proposed by Joseph [6] and Kubota et al. [7] to describe long internal 
gravity waves in a stratified fluid with finite depth. It is written in the form 

ut + -^u x + 2uu x + T(u xx ) = 0, (1) 
where T(-) is the singular integral operator given by 
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T(u)=P I -cot ^{Z-x) u(0dZ (2) 



(P represents the principal value of the integral). Depending on the parameter S (which controls 
the depth of the internal wave layer) (1) reduces to the Korteweg-de Vries (KdV) equation as S — > 
or to the Benjamin-Ono (BO) equation as 5 — ► oo. 

In this paper, we present a new hierarchy which includes the ILW equation at the lowest order. 
This hierarchy is thought to be a novel reduction of the 1st modified KP type hierarchy. In section 
2, we propose a set of differential-difference dispersion relations, and introduce a corresponding 
1st modified KP type hierarchy. In section 3, this 2 + ldimensional hierarchy will be reduced to a 
1 + ldimensional hierarchy which contains the ILW equation at its lowest order. In section 4, we 
discuss more general dispersion relations of differential-difference type. 
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2 The dispersion relations and their corresponding hierar- 
chy 

Let us introduce a pseudo-differential operator, 

w = i- ^d- 1 + w 2 d- 2 + w 3 d- 3 + w 4 d- 4 + ■■■ + w rn d- m , (3) 

where U and Wj (j = 1, 2, ■ ■ ■) are functions of continuous variables t = (ti,t 2 , ■ ■ •) and a discrete 
variable z. We sometimes use x instead of t\ to follow the convention. d~ n denotes 

*~ (4) 

If we employ the Leibniz rule, 

= £ ^ d n-r, (5) 

then d n can be a well-defined operator even for negative integer n. Though the theory is developed 
for the case of m — > oo in general, we, in this paper, confine ourselves to (3) for simplicity [2]. It 
is remarked that the essence of the general theory is still kept in this simplification. 
Let us consider the ordinary differential equation, 

Wd m f(t, z) = (d m - ^S™" 1 + w 2 d m - 2 + ■■■ + w m )f(t, z)=Q (6) 

which has m linearly independent solutions, f^(t,z),f^ 2 \t,z),---,f^ m \t,z). We assume here 
that f^\t, z),f( 2 \t, z), ■ ■ ■ , f( m \t, z) satisfy the following dispersion relations, 

•l-rifl f(j) on f(j) 

A z /Si d fUi CJ = 1, 2, • • • , m and n = 1, 2, • • •), (7) 

where A z denotes a difference operator, 

g(z + 2iS)-g(z) W 
2iS 

(S is an arbitrary constant). 

It should be remarked that U and Wj are expressible by means of a r-function, which is the 
Wronskian of z),f 2 \t, z), • • • , f(- m \t, z), i.e. 

TT n T x 1 (T~t 2 T xx \ 

U = 2-,w 2 = -{--—) (9) 
Let B n (n = 1, 2, • • •) and C be pseudo-differential operators, 

B n = (wd n w- 1 )+, (10) 
C = (WdW- 1 )+, (11) 

where W(z) = W(z + 2i5) and (A) + denotes the differential part of the pseudo differential operator 
A. We use 7 to denote a shift operator: z — > z + 2iS . Then we introduce time evolution equations 
for W by 

, dW 

i — = B n W — Wd n 

W W - (- = 1 ' 2 '-") J (12) 

— = CW - Wd 

2io 
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which are Sato type equations. From (12) we get an infinite system of the Zakharov-Shabat type 
equations 



\C)t 



Bi-Bi 
2i,S 



CB, - B,C = 



i 1 - l {B k ) tl -i 1 - k {B l ) tk + \B k ,B l ]=0 

(M = l,2,---). 



(13) 



Furthermore from (13) we can deduce a system of partial differential-difference equations for U. 
The first few are explicitly given by 



i(U - U) t2 + -(fj-U) x + (U- U)(U - U) x + {U + U) x 

3? - 1 - 

1{U - U) t3 - -{U - U) xt2 + -{U- U) xxx 



0. 



Zi 



- -(f- r ).,-.,■ + ^(u-u)(u-u) x --^s(u-u) x 

3,- „ s „ N Z, r - r _ ,- 3, 
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+ ^(U + U) XX (U -U) + -{U + U) X {U - U) x + -{U - U) 2 (U - U) x = 0, 

-i{U - U) u -\{u + U) t2t2 + l -{u - U) xxt2 + {U + U) xxxx 

+ Z(U - U) X {U - U) xx - l -(u + U) t2 (U - U) x - i{U - U) t2 (U + U) x 

-\{U + U) xt2 {U -U) + ±(U- U) t2 + ±(U + U) xt2 + l -{0 - U) xxx 

+ -(U-U)(U- U) t2 + {U-U)(U- U) xxx - l -{ti - U) 2 (U - U) t2 = 0. 
o 2 

Substituting (9) into (14)-(16), we also have the equations for r by 



iD t2 + -D x + Dij?-T = 0, 

Zi 3 

4D t3 - ZiD x D t2 +Dl + jD 2 x - ^D x ) f • r = 0, 



-2iD 



Dl - iD t2 D 2 x 



1 i 
-^D t2 D x + pD t2 | 



(14) 



(15) 



(16) 



(17) 
(18) 
(19) 



(D denotes Hirota's differential operator). It should be remarked that (17)-(19) are essentially the 
same as the first few of the 1st modified KP hierarchy [8], 



(£>£ +D t2 )r n -T n+1 = 0, 

{Dl - AD t3 - ZD tl D t2 )T n • r„+i = 0, 

{-D 2 t D t2 + D 2 t2 + 2D u )r n ■ r n+1 = 0. 



(20) 
(21) 
(22) 



3 Solutions and special reductions 

iV-soliton solutions for (17)-(19) will be written in the form 



1 + c 1 e v( - t ' Pl ^~ ri ( t,qi ^ 
Kqi) N ~ X + l{pi) N - 1 c 1 e r >( t 'Piy- r K t >n) 



1 + CNe n( - t ' PN ^ +v< - t - qN ' > 
Kin) + l{pN)c N e Tl(t ^ N) ~ Tl(MN) 
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= z n< 

Jci \jeJ 



n 



a jj' 



\3<3' / 



exp ( ^2v{t,Pj) - v(t,Qj) 



(23) 



where the summation is taken over all subsets J of I = {1,2,- • - ,m}. r](t,p), l(p) and 
defined by 



n=l 
„2i<5p _ i 

" -5B-- 

_ %) - /(pjQ %) - %/) 



KPj) - Kir) Kij)- l (Pj') ' 

and the Cj are constants. Notably, the bilinear identity to r is given as follows. 
LEMMA 3.1 For arbitrary t = (ti, t<i, • • ■)> f = (t[,t' 2 , ■ ■ •) and z, r satisfies 

I e -i€(fc,t-t')( 1 + 28k) T {z + 2i6, t - ie(fc- 1 ))r(z, t + ^(fc- 1 ))^- = 0, 
J 2iri 



(24) 

(25) 
(26) 



(27) 



where £(k,t) = X)^=o^™^™ > = (i> 2p"> ' ' ' ' ^F"' ' ' ')• ^ e curve is fafcen around oo and 

excludes the singular points il(p),il(q). 

Proof. Substitute (23) into (27) and we see that Res(il(p n j) + Res(il(q n )) = for Vn. 
We impose the condition 



p j -q j = l{p j )-l{q j )=k j 



□ 

(28) 



on (23). Then, it is reduced to 





( 


\ 


(s«) 








\ 3<3' 


J 



exp ( kjz + n n (kj)t r , 



(29) 



71 = 1 



for appropriately-defined functions /J, n (k). It should be noticed that (29) is the same as the soliton 
solution of the ILW equation [9]. We here present the 1-soliton solution as an example. 



r = 1 + c exp 



kz + kti + (k 2 cot kS - -)t 2 



4 ( k 3 - 3fc 3 cot 2 kS + ^k 2 cot kS - -Jrfc ) t 3 

' 1 Q z J 

3 13 1 \ 

-k 4 cot kS + k 4 cot 3 kS - -k 3 cot 2 kS + -k 3 + —k 2 cot kS — —k ) t± H 

o z A I 



(30) 



We can regard z and t\ as the same variable under this reduction, because the coefficient of z is 
equal to that of t\ in the exponentiated part of (29). It also should be noticed that 

LEMMA 3.2 If kj 's are real, ^ n (kj) and ajj' are also real (j,n= 1, 2, • • •). 

Proof. 

From (25) and (28), we get 



l(pj) + l(qj) = —ikj cot kjS + 



(31) 



Because l(pj) + l(qj) is purely imaginary and l(pj) — l(Qj){= kj) is real, there exist real r, 9 by 
which 



l(p,)=re M , l(q 3 )=re^- e \ 



(32) 
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By the definition of n n {kj), we have 

H n (kj) = i"- 1 (l( Pj ) n - i( qj ) n ) 



i 



{ 



-1 (r e in6 _ re in(n-e)^ 

2i n ~ 1 r n cos(n9) for n odd 
2i n r n sin(n#) for n even. 



(33) 



Hence fi n (kj) is real. We also deduce 



Kpj) - Kir) Ku) - l (Pj') 

{re w ~r'e ie '){re- %e -r'e- 16 ') 

r 2 + r /2 _ rr / cos (0_0') 



(34) 



r 2 + r' 2 — rr' cos(0 + 6') 

which gives that ajji is real. rj 
By this reduction, the x-shifts can take the place of the z-shifts and (14)-(16) are rewritten 
into the equations for U (x) := U(x =p iS), i.e. 

i(U~ - U+) t2 + l -(U~ - U+) x + (U- - U+)(U- - U+) x + (U- + U+) xx = 0, (35) 
2(17" - U + ) u - ^(U~ - U+) xt2 + hp- - U+) xxx + ^(U~ + U+) xx (36) 

+ j(u~- u+)(u- - u+) x -^(u~- U + ) x + |(ET + u+) xx {u- - U+) 
+ l(u~ + u+uu- - u+) x + \{u~- u+f(u- - U+) x = 0, 

-i(U~ - U+) u - 1(17- + U+) t2t2 + l -(U~ - U+) xxt2 + (U- + U+) xxxx (37) 

+ 3(tr - u+uu- - u+) xx - l -(u~ + u+) t2 (u- - u+) x - i(u~ - u+) t2 (u- + u+) x 

- \{U~ + U+) xt2 (U- - U+) + ±(U~ - U+) t2 + h {U - + U+)xt2 + 1 _ {U - - U+) xxx 

+ -AU-- u+)(u- - u+) t2 + (u- - u+){u- - u+) xxx - Uu- - u+f{u- - U+) t2 = 0. 

i 

If we consider lemma 3.2 and suppose that U(x) is analytic in the horizontal strip between 
Imi = — i5 and Imi = i5, we can introduce a dependent variable u which satisfies [9] 

u= l -(U--U+), (38) 

T{u) = \{U- + U + ). (39) 
Substituting this u into (35)-(37), we obtain 

u t2 + \u x + 2uu x + T(u xx ) = 0, (40) 
o 

-4it t3 - 3T(u xt2 ) - u xxx - 6uT(u xx ) - 6u x T(u x ) + I2u 2 u x (41) 

12 3 3 

— — uu x + -T{u xx ) + -^u x = 0, 

-2u u - T(u t2t2 ) + u xxt2 - Au 2 u t2 - 2u x T{u t2 ) - Au t2 T{u x ) (42) 
2uT(u x t 2 ) + 2T(u xxxx ) — \2u x u xx ^uu xxx 

x 2 4 1 

+ -^T(u xt2 j + -u xxx - -uu t2 + -pU t2 = 0. 
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Because the lowest order is nothing but the intermediate long wave equation, this hierarchy should 
be called the ILW hierarchy. 



4 More general dispersion relations 

Other possible differential-difference dispersion relations than (7) can be 

f il-n^jU) =d n f (j) 

1 i s ~ l A z f^ =0V (J '\ 



(43) 



where s is some fixed positive integer, j = 1, 2, • • • , m and n — 1, 2, • • •. Notably, (43) corresponds 
to (7) when s = 1. 

For each s, we can deduce the corresponding hierarchy in the same way as in the preceding 
sections. Hence we just present the bilinear identity which generates the hierarchy for r [8], 



I e -«(fc,t-f )(! _ 25(-k) s ) T {z + 2i6, t - ^e(fc- 1 ))T(z, t + le^ 1 ))^- = 0. 

J 27TI 



(44) 



The ./V-soliton solution to the hierarchy is written in the form 



1 + c 1 e H( - t < Pl ^ H ( t > qi '> 
L( qi )+ L(p 1 )c 1 e fl ( i *)- ff (^') 

Hq^- 1 + L{p 1 ) N - 1 c 1 e H ( t 'P^- H< - t -'^ 

/H(L( qj ,)-L( qj )) 

j'>j 

e (n c i] n A i? I cx p ( e - <&) 



1 + c N e H{t - PN)+H{t > qN) 
L{q N ) + L{p N )c N e H ^-P^" H ^^ 

L(qpf) N ^ 1 + L{p N ) N - 1 c N e H ^ t -P^' H{t ^ 



(45) 



where the summation is taken over all subsets J of I = {1, 2, • • • , m}. H(t,p),L(p) and A,y are 
defined by 



ff(t,p)=pz + ]T*"- 1 L(p)'X, 

n=l 
e 2 l 5p _ 1 



i l - s L{pf = 



2iS 



(46) 

(47) 
(48) 



" L( Pj )-L(q r ) L( qj )-L( Pj ,) ■ 

It should be noticed that L(p) is multi- valued. It is easy to check that this soliton solution satisfies 
the biliniear identity (44). 

We impose the reduction condition, 



Pi-oj = i 8 ' 1 {L{ Pj ) s - L( qj y) 
= (-ly- 1 - i — 

= kj (j = l,2,...). 



Substituting (49) into (45), we have 

t = E (li C i) II A n' CX P I E ( k > z + E M n(kj)t n 



(49) 



(50) 



Jd \je.J 



\3<3' ) 



for appropriately-defined functions M n {k). We can regard z and t s as the same variable under this 
reduction. Furthermore, the following lemma holds. 
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LEMMA 4.1 If kj 's are real, M n (kj) and Ajji are also real (i = 1, 2, • • ■). 
Proof. ^From (48) and (50), we see that 

L( Pj y - L( qj y = (if-'kj, (5i) 



(52) 



If s is odd, L(pj) s — L(qj) s is real and L(pj) s + L(qj) s is purely imaginary. Hence there exists real 
r, < 9 < 2tt by which 

L{p,) s = re 16 , L(q 3 ) s = re l ^- ff) . (53) 
If s is even, L(pj) s — L(qj) s is purely imaginary and L( P j) s + L(qj) s is real. Then we have 

L( Pj y=re ie , L( qj y=re- ie . (54) 
For both cases, we choose L(pj), L(qj) as 

L( Pj ) = ^, L(q 3 ) = ^e^-^. (55) 

Because L(pj) — L(qj) is defined to be real and L(pj) + L(qj) to be purely imagenary, we see that 
Mj,Ajf are real by means of lemma 3.2. This, at the same time, takes care of the multi-valuedness 
of the function L(p). rj 
If we consider lemma 4.1 and suppose that U(t s ) is analytic in the horizontal strip between 
Im t s — —i8 and Im t s — i5, we can introduce a dependent variable u which satisfies 

u= l -{U{t s +i5)-U(t s -i5)), (56) 

T s (u) = ^(U(t s +iS) + U(t s -iS)), (57) 

where 

/°° 1 

Ys cot &^ ts) \ u{0dt (58) 

As another concrete example, different from the ILW hierarchy, we apply the preceding argu- 
ment to the case s = 2. From (44), we get 

I e -«(fc,t-0 n _ 25k 2 ) T {z + 2^<S, t - ^e(fc- 1 ))T(z, t + J£ (fc- 1 ))— = 0, (59) 
J 2m 

which generates the hierarchy for r. The first few equations are 

-2D t3 +D 3 X + 3iD x D t2 + ^D x ^j f ■ r = 0, (60) 

(&D U - AiD t3 D x - ZiDl - %D% - ^D 2 X - | D t )j f ■ r = 0. (61) 

This hierarchy is essentially the same as the 2nd modified KP hierarchy [8] 

(2D t3 +Dl + 3D x D t2 ) r„ • r„ +2 = 0, (62) 
(Df - 4 J D 1 L>3 - 3,D 2 2 - 6£> 4 ) r„ • r„ +2 = 0, (63) 

The 1-soliton solution for this hierarchy is written in the form 

T = l + Cexp[H(t,p)-H(t,q)], (64) 

OO 

H(t,p) = (p-q)z + J2^ 1 L(p) n tn, (65) 



n=l 



2 -e M P + 1 



W = (66) 
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By the reduction condition 

p-q = l{L(p) 2 -L(q) 2 ) 
e 2iSp _ e 2i6q 



2iS 

= k, (67) 



(64) is reduced to 



t = 1 + cxp 



kz+Jk cot kS - S- 1 + cot kS - 5- 1 ) 2 + k 2 t 1 + kt 2 (68) 



+ \jkcotk5 - S- 1 + VWcot kS - 5- 1 ) 2 + k 2 

x (-fccot kS + S- 1 + ^/(k cot kS - S- 1 ) 2 + k 2 /V2^j t 3 + ■ 

Thus the t 2 -shifts take the place of the z-shifts. If we impose the analytic condition as before on 
U, there exists u which satisfies 

u= l -{U{t 2 +i5)-U{t 2 -i8)), (69) 
T 2 (u) = ^(U(t 2 +iS) + U(t 2 -i6)). (70) 

Now, (60) and (61) are expressible by means of u as 

2u t3 - u xxx + 3T 2 (u xt2 ) - 3 (uT 2 (u x )) x + 3u 2 u x (71) 
— 3[u[ Ut 2 dx ) — —u x = 0, 



S 

6u t4 + T 2 (u xxxx ) - A{uu x ) x + 12T 2 (u)T 2 (u) x + 4u 3 u x (72) 
-6 (6u 2 T 2 (^)) x + AT 2 {u xt3 ) - Au x f u t3 d£ - Auu t3 + 3T 2 (u t2t2 ) 

J — OO 

-6u t2 J u t2 d£ + —T 2 (u xx ) - —uu x - -u t2 = 0. 
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